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Abstract. We look at multiparticle production processes from the Information 
Theory point of view, both in its extensive and nonextensive versions. Exam¬ 
ples of both symmetric (like pp or AA) and asymmetric (like pA) collisions are 
considered showing that some ways of description of experimental data used in 
the literature are of more general validity than usually anticipated. 
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Introduction 

The multiparticle production processes are usually hrst approached by means of 
statistical models [ ^ in order to make quick estimations of such parameters as 
temperature T or chemical potential p of the hadronizing matter (with tacit as¬ 
sumptions that T and p have the usual meaning in the realm of hadronic produc¬ 
tion). It means that ’’thermal-like” (i.e. exponential) form of relevant distribution 
(in transverse momentum pT or in rapidity y) is used to fit data, 

/(A)^exp[-A/A]. (I) 

This formula occurred to be very robust, mainly because (cf., [ 12] ) the N — 1 
unmeasured particles act as a heath bath which action on the observed particle is 
described by a single parameter A identihed (for system in thermal equilibrium) with 
temperature T. Discrepancies from (D observed in many places are then attributed 
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to the fact that in reality such a ’’heath bath” is neither infinite nor homogeneous 
and one needs more parameters for its description. The minimal extension would 
be to regard A as being X-dependent, for example, A ~ Aq + a • X [|2|, what results 
in 

/g(X) ~ exp^ [-X/Ao][1 - (1 - where 9 = 1 + a. ( 2 ) 

In the case when X is energy E and A is temperature T the coefficient a is just 
the inverse heat capacity, a = IjCv [El- Such phenomenon leads to the so-called 
nonextensivity with q being the nonextensivity parameter. Notice that for (7 —> 1 (or 
for a ^ 0) fq{X) becomes /(X). We show here that eqs. ([Q and m originate in 
a natural way from information theory (IT) approach in its, respectively, extensive 
(based on Shannon entropy) and nonextensive (based on Tsallis entropy) forms. 
In what follows (Section J we shall briefly present this approach illustrating it (in 
Sections □ and p) with some fits to the existing data. 

Information Theory and Multiparticle Production Processes 

To introduce IT in the present context let us consider typical situation: experi¬ 
mentalists obtain some new and intriguing data. Immediately these data become 
subject of interest to theoreticians and in no time a number of distinctive and unique 
(as concerns assumptions and physical concepts) explanations is presented, which, 
disagreeing about physical concepts used, nevertheless all fit these data. The nat¬ 
ural question arises: which of the proposed models is the right one? The answer 
is: all of them (to some extent). This is because experimental data are providing 
only limited amount of information and all models mentioned here are able to re¬ 
produce it. To quantify this reasoning one has to define, using /T“, the notion of 
information. Its extensive version is based on the Shannon information entropy, 

Inp,, (3) 

i 

where pi denotes probability distribution of interest. The least possible informa¬ 
tion, corresponding to equal probability distribution of N states, pi = 1/N, results 
in maximal entropy, S = InX. The opposite situation of maximal information, 
when only one state is relevant (i.e. pi = 1 and = 0 ) results in minimal en¬ 
tropy, S = 0. Denoting by {Rk) « priori information available on experiment, like 
conservation laws and results of measurements of some quantities Rk, one is thus 
seeking probability distribution {pi} such that: i) it tells us the truth, the whole 
truth about our experiment, i.e. in addition to being normalized it reproduces the 
known results: 

^Pi = I and Pi Rk{xi) = {Rk) , (4) 

i i 

and ii) it tells us nothing but the truth about our experiment, i.e. it conveys the 
least information (only the information contained in this experiment). 
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To find such {pi} one has to maximize the Shannon entropy under the above 
constraints (therefore this approach is also known as MaxEnt method). The resul¬ 
tant distribution has familiar exponential shape 



^ ^ A/e ■ 
k 


( 5 ) 


Although it looks identical to the ’’thermal-like” (Boltzmann-Gibbs) formula (^3) 
there are no free parameters here because Z is just normalization constant assuring 
that A/e are the corresponding lagrange multipliers to be calculated 

from the constraint equations 0.^ 

It is worth to mention at this point [ 0] that the most probable multiplicity 
distribution P(n) in the case when we know only the mean multiplicity (n) of 
distinguishable particles is geometrical distribution P{n) = (n)”/(l -1- 
(which is broad in the sense that its dispersion is D{n) ~ (n)). Additional knowledge 
that all these particles are indistinguishable converts the above P{n) into Poissonian 
form, P{n) = (n)”exp(—(n))/n! which is the narrow one in the sense that now its 
dispersion is D(n) ^ \/(^. In between is situation in which we know that particles 
are grouped in k equally strongly emitting sources, in which case one gets Negative 
Binomial distribution [[?]'’’ 


P[n) = r{k + n)/ [r(n -k i)r(fc)] 


/ k 


k ' 


/ 

1 + 7—r 

\{n). 

L (’^)J 


The other noticeable example provided in [|21 is the use of IT to find the minimal 
set of assumptions needed to explain all multiparticle production data of that time. 
They were equally well described by models like multi-Regge, uncorrelated jet, 
thermodynamical, hydrodynamical etc., which, after closer scrutiny, turned out to 
share (in explicit or implicit manner), two basic assumptions: i) that only part 
W = K ■ s/s (0 < if < 1) of the initially allowed energy is used for production 
of observed secondaries (located mostly in the center part of the phase space; in 
this way inelasticity K found its justification [ini, it turns out that K ~ 0.5); ii) 
that transverse momenta of produced particles are limited and the resulting phase 
space is effectively one-dimensional (dominance of the longitudinal phase space). 
All other assumptions specific for a given model turned out to be spurious. 

Suppose now that some new data occur which disagree with the previously 
established form of {pi}- In IT approach it simply signals that there is some addi¬ 
tional information not yet accounted for. This can be done either by adding a new 
constraint (resulting in new A/^+i, cf., for example [|^, we shall not discuss it here) 
or by using some other form of /T, for example its nonextensive version (ITq). The 
later is necessary step for systems which experience long range correlations, mem¬ 
ory effects, which phase space has fractal character or which exhibit some intrinsic 
dynamical fluctuations of the otherwise intensive parameters (making them exten¬ 
sive ones, like T here). Such systems are truly small because the range of changes 
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is of the order of their size. In this case the Shannon entropy © is no longer a 
good measure of information and should be replaced by some other measure. Out 
of infinitely many possibilities [cni we shall choose Tsallis entropy m, 

-S'? = ^ Sq=l=-J2p^\np„ (6) 


which goes over to Shannon form © for (7 ^ 1. The {p\'^^} are obtained in the 
same way as before but with modified constraint equation:® 

E Rk{x.) = (7) 


This leads to formally the same formula for pi = p^^^ as in Eq. © but with 
Z ^ Zq and exp(...) ^ exp^(...). Because such entropy is nonextensive, i.e., 
Rq{A+B) ~ RqA T SqB T (1 Q^^qA ' SqB [[OJ the whole approach became known as 
nonextensive (Tsallis) statistics. It should be stressed at this point that the nonex- 
tensivity parameter q cumulates action of all possible dynamical sources causing 
deviation from the usual Boltzmann-Gibbs statistics or Shannon entropy and as 
such can be considered as a useful phenomenological parameter allowing to de¬ 
scribe data without deciding which of dynamical models leading to the same q is 
the right one [Cl- However, in what follows we shall mainly address one possible 
sources of g ^ 1, namely intrinsic fluctuations present in the system represented by 
fluctuations of the 1/A parameter. This is because, as shown in [E]j in this case 


g = 1 ± 



(1/A)^ /(1/A)^ 


( 8 ) 


i.e. parameter g is a measure of such fluctuations (with (...) denoting the respective 
averages over them.7) 


Confrontation with Experimental Data - Symmetric Case 

We shall now confront these ideas with reality. At first we shall remind shortly 
main results of our recent investigations of hadronizations taking place in collisions 
of symmetric systems like pp and pp [II211IS11TH1, heavy ions AA 113 cm EH or e+e 
[E2- This will be followed by some new results on collisions of asymmetric systems 
exemplified by pA collisions. It must be stressed that what we are proposing is not a 
new model but rather sets of least biased, most probable distributions describing data 
by accounting for available information provided in terms of constraints emerging 
from conservation laws and from some previously known experimental facts (or from 
some assumed dynamical input, which is thus confronted with experimental data). 

Hadronization means that some invariant energy M (assumed to be known) 
gives rise to a number N of secondaries (also assumed to be known) and question 
asked is: in what way these secondaries are distributed in the allowed phase space? 
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So far distributions in transverse momenta pt were treated separately [El ED from 
distributions in rapidity [El ^1 El El 1^ - for which it was always assumed that 
mean transverse mass px = -\/m^ + (pt)^ was known, i.e. that we are working in 
the effectively 1-dimensional phase space. 

Concerning px, distributions it has been shown that using q-statistics (with 
g ~ 1.05 for AA collisions [ED to g ~ 1.1 for pp collisions [El) one can describe data 
in large domain of px, see, for example, Fig. ^ The characteristic feature is that 
now values of q are much smaller than those obtained fitting data in longitudinal 
phase space (cf. [El for discussion). Following [El[^ (see also [El El El) we argue 
that q > 1 shows existence of the fluctuations of temperature in the hadronizing 
system mentioned at the beginning. In fact, the q = 1.015 for AA collisions [EH 
corresponds (according to ©) to fluctuation of T of the order AT/T ~ 0.12, which 
do not vanish with increasing multiplicity [El- These fluctuations exist in small 
parts of the hadronic system with respect to the whole system, they are not of the 
event-by-event type, for which AT/T ~ 0.06/^/N —s- 0 for large N. It should be 
stressed at this point that such fluctuations are very interesting [ESI because they 
provide a direct measure of the total heat capacity of the system, C: 


g^(/3) dg ^ _ 

W c 


(9) 


In fact, because C grows with the collision volume V of reaction we expect that 
( 7 (hadronic) > ^(nuclear) which seems to be indeed observed (cf. also [^J where 
nuclear collisions data at different centralities providing direct access to volume V 
were analyzed). 

When applied to rapidity distributions IT method leads to formula formally 
identical with formulas used in statistical models (cf., for example, like [I24|b 

^ pyM 

fN{v) = ^ • exp [—/3 • cosh y] where Z= I dy exp [—/3 • cosh y]. (10) 

^ J-Ym 

However, whereas in [E] (and in other similar models) 1/Z and jS were just two 
free parameters, here Z is normalization constant and /3 = f3{M^ N, px) is obtained 
by solving constraint equation,® 

/ Ym pYm 

fniy) = 1 and / dy [px ■ coshy] • /Ar(y) = —. (11) 

-Ym J-Ym 


It means that the parameter (3 (inverse of the so-called partition temperature intro¬ 
duced in [1^ type of models) is connected (via IT method used) with the dynamical 
input given by: the allowed energy M (usually taken as a fraction, M = K^/s^ of 
the total energy of reaction with parameter K G (0,1) being the so called inelas¬ 
ticity of the reaction), number of secondaries produced N and the mean transverse 
mass Px- In asymmetric collisions to be discussed in the next section one would 
have to add also constraint imposed by momentum conservation (which is satisfied 
automatically in the case of symmetric collisions discussed at the moment). 
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As detailed in [H ns, fN{y) changes from fN= 2 {y) = ^[S{y-YM)+S{y + YM)] 
(with P(N —> 2) ^ — oo) via /nq = const, (with f3{No ~ lii{M/fj,T)) — 0) to 
fN^N^ax = ^(y) (with l3{Nmax = M/^t) = +oo). In other words, for small 
multiplicities /3 is negative (a feature alien to any statistical model!) and it becomes 
zero only for N Nq. At this multiplicity fN{y) = const., behaviour known 
as Feynman sealing.^ Its occurence means that energy dependence of multiplicity 
follows that of the longitudinal phase space. For N > Nq additional particles have 
to be located near the middle of phase space in order to minimize energy cost of 
their production. As result /3 > 0 now, in fact (see [0] for details) for some ranges 
of {E) = M/N quantity P = P{E) remains approximately constant. Needless to 
say that for N N^ax oil partieles have to stay as much as possible at the center 
and therefore /3(M, N N^ax) +oo, whereas fN^Nmaxiy) ^{v)- 

For nonextensive approach exp(...) must be replaced by expg(...)q and Nq 

~ (2\nNmaxY■ As shown in f I12| q acts now as a free parameter allowing 
for changing separately the shape and the height of (which were interlocked 
for q = 1 case). For q > 1 one enhances tails of distribution reducing at the same 
time its height. For g < I the effect is opposite and additionally there is kinematical 
condition, 1 — (1 — q)PqfiT coshy > 0 , reducing in this case the allowed phase space.® 




P, [GeV] y 

Fig. 1. Example of description of data on pT spectra from UAl experiment on pp collisions [ 
13 (left panel) and data on rapidity spectra in pp and pp collisions \ U2j (right panel). 

In [E] we have used to fit pp and pp data introducing inelasticity K as 
a free parameter to be deduced from data (cf. Fig. P). We shall not discuss here 
the inelasticity issue (see [ [HI for details and further references) but concentrate 
on the q parameter, which turned out to follow the same energy dependence as the 
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experimentally deduced parameter Xjk of the Negative Binomial (NB) multiplicity 
distribution [ESI- This finding prompted us to argue that parameter q in this 
case is accounting for a new bit of information so far unaccounted for, namely for 
fluctuations in the multiplicity of produced secondaries (notice that in the constraint 
equation for (3 we were using in this case experimental values of the corresponding 
mean multiplicities only). In fact, it is known [Ej that NB can be obtained from 
Poisson distribution once one allows for fluctuations in its mean value n of the 
gamma distribution type, namely 

r(fc + n) 7^ .,2i 

r(l + n)r(A:) ■ (7 + 1)'=+" ^ ’ 

^ _ cr'^inch) _ 1 

k {uchy {rich)' 


P{n) = 


dn 




m 


where 


7 = 


and 




Fig. 2. Example of description of data on rapidity spectra from NA49 (left panel) for negatively 
charged pions produced in central PbPb collisions at different energies [13; the best fit for 17.3 
GeV is with additional (assumed) information of existence of two extensive sources located at 
y = diAy = 0.83 in rapidity (see text for details). Right panel: fits to PHOBOS data for the most 
central Au + Au collisions [HSl 

Assuming now that these fluctuations contribute to nonextensivity defined by 
the parameter q, i.e. that D{n) = q — 1 one gets 

q = I + —, (13) 


what we do observe (cf. also [|H|). 
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In Fig. [3 we show comparison with AA data at SPS (NA49) [|S| and RHIC 
(PHOBOS) [EOI energies. The NA49 data can be fitted with q — 1.2, 1.16 and 1.04 
going from top to bottom and so far we do not have plausible explanation for these 
number. However, it should be stressed that at highest energy two-component 
extensive source is preferable. The PHOBOS data with q = 1.27, 1.26 and 1.29 
going from top to bottom. Extensive fits do not work here at all. The first clear 
discepancy has been found when fitting e+e“ data (see [ 122 ]) where dn/dy with 
clear minimum for y = 0 cannot be reproduced in such simple ITq approach. It is 
thus obvious that there is some new information we did not accounted for. It looks 
like we have two sources here separated in rapidity (two jets of QCD analysis) but 
of no statistical origin (rather connected with cascading process) [| 22 |- 

Confrontation with Experimental Data - Asymmetric Case 

Let us now proceed to a more complicated case of asymmetric collisions exemplified 
by pA processes [l2nil2ZII2HI In this case: i) both colliding objects are different and 
have different masses (one can therefore expect that energy transfer to the central 
region from each of them is not necessary the same as it was assumed before); ii) 
the pA collisions introduce a new element of uncertainty, the effective size of the 
target, i.e. the number of nucleons, v, from the nucleus A with which the impinging 
nucleon is interacting and the way this interaction proceeds.-’ The number v can 
be estimated from measurement of the number of ’’gray tracks” [ 1221 , therefore it 
will be included to our input information. We shall analyze data [ 122 ! in which 
attempt was made to isolate dN^/dy for v = 1, 2, 3,4 and data [[22] where only the 
mean v number of collision is known together with dNp/dy. The a priori available 
invariant energy in such a case is equal to (we neglect all nuclear binding and 
Fermi motion effects): = vs + {v — l)^m^ where s = 2m? + -I- is 

the invariant energy squared for NN collisions. 

The assumed knowledge of v must be supplemented by the known dependence 
of total mean multiplicity of secondaries produced in the central region, on the 
mean number of inelastic collisions {v ), [1221 EZj: 

K = lil + {v))N (14) 

{N is multiplicity of particles produced in collisions at the same energy y^s,^=i)A 
To be able to apply IT methods we must additionally decide whether pA collision is 
more like a two body collision between a kind of ’’tube” of mass mi, = vm contain¬ 
ing V collectively acting nucleons and single nucleon of mass m (cf. Fig. |2Ia) and 
(b)), or whether it is rather a sequence of v consecutive collisions of the impinging 
nucleon with v nuclear nucleons (cf. Fig. I^:). 

In the first case we have situation similar to considered before with the following 
formula for the rapidity distribution of secondaries produced in two-body p{vN) (in 
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(b) 


(E„-P,) 

(E®-pf) ce;,-p;)^ 

(m 

(Ei,Pi) 1 

(ei, pi) 1 

) • • • 

(Wi, Vi) 1 

(ei, pi) 1 

M ■ ■ (m 

(w,,Vi) 1 

(e., p.) 

J) 

(w„, V,) 


1 


: ^ 


(E., P.) 

(E', p;) 


(c) 


Fig. 3. Schematic views of tube (a) and sequential (c) models of the pA collisions. The case (6) 
is essentially a particular realization of the tube type of model (a). 


the CM frame of the NN system): 


dN, 

dy 

Z. 


My) 



= — • exp sinh y — cosh y ), 

Zv 

dy exp (—Aiz/tt sinhy — (d^yT coshy). 


(15) 

(16) 


We have now two lagrange multipliers, and which are given by the corre¬ 
sponding energy and momentum conservation constraints: 



— pT [Ai/ sinh y—(3t, cosh y\ 


coshy 
sinh y 


Z, I W, = (iyE + K)^, 1 

X^i^yr \p^ = -{vR - K) J ■ 


(17) 

The energy transfer from the projectile nucleon, characterized by inelasticity if, 
is allowed to differ from energy transfers from the nuclear nucleons (cf. Fig. I^b) 
characterized by inelasticities Ri (for simplicity we shall assume Ri = R in what 
follows). The invariant mass Mi, of hadronizing system is now equal to 


- Pi = ^ KRvs P {vR-Kfw? 


R=K 


Ky/s^ (18) 
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whereas the (longitudinal) phase space in which particles are produced is given by: 

^max — — <5jy and = —Yvm — 5^, (19) 


{Yyrn is the same as Ym before but calculated for M^) where 6u being the rapidity 
shift between the NN and {vN) — N center of mass frames: 


tanh 6^ 







( 20 ) 


The apparently asymmetric form of rapidity distribution as given by Eq. 03 is, 
however, an artifact connected with our choice of the reference frame. Changing 
variables in eqs. (iisj-ini from y io y = y -\- i.e. proceeding from the CMS of 
NN to the CMS of {yN) — N, one gets similar distribution as before (i.e. depending 
only on one lagrange multiplier j3v): 

My) ^ My = y-s,) = ^- coshy. ^ (21) 

pYum ^ y 

/ dycoshye-^'^'^-^™"*'^ = ■ M^, (22) 

N„PLT 

where = Mi, {K, R). It means that in such an approach we always can find 
frame in which rapidity distributions f^, are symmetric function given by (I21II . In 
Fig. H we show examples of fits to some available data (notice that data display 
clear asymmetric character which cannot be reproduced by the method used here). 


In the second case pA scattering is assumed to proceed via sequence of the 
1 / two-body processes, cf. Fig. Oh. The resultant rapidity distribution f„{y) is 
then composed of ’’elementary” distribution functions describing collisions of the 
impinging nucleon with the subsequent nucleons in the target nucleus: 

V 

fM) = ^h{y)- (23) 

Contrary to this formula stresses not the fact that in the pA reaction one has 
v+l participating nucleons but that one has here v consecutive collisions treated as 
’’elementary” ones. This can be visualized rewriting (O as N„ = Ni + {v- 1)A^2/2 
where Ni^2 are such that 2iVi -1- {v—l)N2 = {i' + l)N. We have thus ly ’’elementary” 
collisions with a possibility that the first can differ from the remain v—l (in general 
they all are different). As result, the elementary inelasticities (i.e. inelasticities in 
subsequent collisions) are not necessarily the same. In this way one approaches as 
near as possible the notion of independent production for which one demands the a 
priori knowledge of energy of each reaction. Also here we shall differentiate between 
the fractions of energies contributed to each fi by the nuclear nucleons, 
and fractions of energies contributed by the impinging nucleon at each collision, 
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y y 

Fig. 4. Example of fits to pp, pAr and pXe data [I23 by means of the collective tube model 
approach as given by Eq. on with R = K (left panel) and with R ^ K for pA collisions (right 
panel). 


^ 1 = 1 ....,!^ In this case the general form of the energy-momenta {ei,pi) and {wi,Vi) 
flowing to the blob Mi in Fig. O: is given by: 


Ci 


Pi 


Wi 


1=1 


-1 1 

1/2 

- 

2 — 1 

n(l-^,)s-4m" 

1 

sD 1 

n(l-X,)s-4m^ 

_l=i 


_i=i 


Vi = f s - 4 to 2 - iy(1 - • s - 4 to2 


(24) 

(25) 


with Mf = [{a + Wi) ; {pi + Vi)]'^. Notice that the real fraction of energy deposited 
by impinging nucleon in, say, second collision , is equal to = K 2 ■ (1 — Ki), in 
general = Ki ■ 0}=! (1 ~ ^l)- We have now that 


My) 



(26) 




'-Yi 


Z^ 

Mpt 


■ M,. (27) 


with Yi calculated in the same way as Ym before with Mi and Ui replacing W and 
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N, respectively. These distribution will be centered (in the CM frame of NN) at 



fC ^ 


i=i 


Ri=Ki=K (i - 1) , 

^ oc—-—In(l-K), 


(28) 


i.e. the expected shift in rapidity is now linear in the number of collisions, not 
logarithmic one as in Eq. CT . However, for our limited values of v this difference 
is practically not visible. The last input needed to apply IT is some estimation of 
multiplicities Ni. This is the most uncertain point here and can be done in many 
model-dependent ways. Here we have simply assumed that Ni = a - Mf (with a and 
c being free parameters) and made use the fact that N = to 

eliminate parameter a and to write Ni = N/ 1 -|- : where c = 0.556 


is obtained from reproducing the observed multiplicity N (actually, logarithmic 
dependence Ni = a + h ■ InM^ would work equally well). In order not to introduce 
too much freedom with the choices of Ki we were proceeding in the following way: 
Ki was fixed by the v = 1 data and then used as then K 2 was fixed by = 2 

data and used as if ,=2 and so on up to K^, which was hxed by = 4 data.* On 
the other hand the fractions of energy deposited by the consecutive target nucleons 
participating in the collision was kept the same for all of them and equal R — R^. 


>, 

o 

z 

a 




y 


y 


Fig. 5. Left panel: Comparison with pA data for u [ESI using tube model (Eq. on and Fig. 
Et) for K ^ R. Notice that data for u = 1 are not compatible with data for pp collisions from [ 
w\ presented in Fig. O Right panel: the same but using sequential model (Eq. and Fig. Et)- 
Ki is the fraction of the actual projectile energy in the i — th collision deposited in central region. 
The fractions of energy deposited by target nucleons is kept the same and equal R = Rj ^. 


Notice that in the sequential model we have clearly smaller consecutive energy 
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transfers (as given by from the nucleon traversing the nucleus (this follows 

observation made in [EOl)- On the other hand, both approaches lead to essentially 
the same results. However, only sequential model is potentially able to reproduce 
the visible asymmetry seen in data for larger v because ” elementary” ji enter with 
different weights and cover slightly different regions of phase space. However, at 
present energies the effect is not as large as expected and sequential model leads to 
essentially the same results as tube model.™ To account for the shift in the position 
of maximum of rapidity distribution in the tube model we have to allow for K ^ 
i.e. for different energy depositions from each of the projectiles. However, even then 
we cannot produce the asymmetry of rapidity distribution seen in data, in particular 
we are not able to correctly describe the part of distribution connected with the 
impinging proton. It should be stressed that using in this case the nonextensive 
approach ITq would not help because also in this case the obtained distributions are 
symmetric in the rest frame of the hadronizing system. There is clearly difference 
between the impinging propton and nuclear parts of the phase space. We conclude 
therefore that, according to the rules of information theory approach, data [ESIEZI 
still bear some additional information, not identified so far. 

Summary and Conclusions 

We have demonstrated that large amount of data on multiparticle production can be 
quite adequately described by using tools from information theory, especially when 
allowing for its nonextensive realization based on Tsallis entropy. We have argued 
that nonextensivity parameter q entering here can, in addition to the temperature 
parameter T of the usual statistical approaches, provide us valuable information on 
dynamical fluctuations present in the hadronizing systems. Such information can 
be very useful when searching for phase transition phenomena, which should be 
accompanied by some specific fluctuations of non-monotonic character [ES|. 

As concerns the asymmetric collisions example of pA data, we have demon¬ 
strated that they contain additional information to the usual one used in Section 
B The immediate candidates are the rescattering of produced secondaries in the 
nucleus (effectively increasing values of and making it y—dependent and in this 
way limiting rapidity space available on the nuclear side) or diffraction dissociation 
part of the production (which can take different shape than in the pp collisions), 
but they can be other possibilities as well [El. 

Finally let us comment shortly on the inelasticities in pA collisions obtained. 
In the tube model they are, in general, increasing (or, at least, non-decreasing) 
with A or with z/, cf. Figs. Em and El Notice that for K ^ R case the effective 
inelasticity, which is of interest for any statistical model approach (i.e. the part 
of the total energy available for production of secondaries in the central region 
of reaction) is K^ff = VK ■ R. In sequential model this inelasticity is clearly 
decreasing with the consecutive collisions (in agreement with what was found in [ 
EOl)- However, the more precise statement could be only done with a data taken 
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at much higher energies. Notice that we did not present fits to pAr and pXe data [ 
EH in this case because we would have to introduce here another piece of additional 
information represented by distribution of number of collisions, which is model 

dependent quantity. The high energy counterpart of data [EH would be therefore 
most welcome. The point is, however, that in addition one should also have data 
taken for pn inelastic collisions (possibly with the help of deuterons) as it is possible 
that at least part of the inconsistency between pp data from [EZl (cf.. Fig. 0 and 
1 / = 1 data from [EH in Fig. Elis because the later contain also contribution from 
pn reactions." 

We close with the statement that results presented here are very encouraging 
and call for further systematic effort to describe existing data in terms of (T, q) for 
different configurations and energies in order to find possible regularities in their 
system and energy dependencies and possible correlations between them. 
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Notes 

a. The complete list of references concerning IT relevant to our discussion and 
providing necessary background can be found in [ElIH- The value of Boltzmann 
constant is set to unity, k = 1. 

b. Notice that using the entropic measure S = J2i [Pi T (1 ± pi) In (1 ± pi)] 
(which, however, has nothing to do with IT) would result instead in Bose- 
Einstein and Fermi-Dirac formulas: pi = {I/Z) ■ [exp[/3(ei — p)] =F 1]~ , where 
[3 and p are obtained solving two constraint equations given, respectively, by 
energy and number of particles conservation [ E] . It must be also stressed 
that the final functional form of pi depends also on the functional form of the 
constraint functions Rk{xi). For example, R(x) oc ln(a;) and ln(l — x) type 
constrains lead to pi oc a;“(l — Xi)^ distributions. 

c. It is straightforward to check that the Shannon entropy decreases from the 
most broad geometrical distribution towards the most narrow Poissonian dis¬ 
tribution. 

d. The most drastic situation was with the multi-Regge model in which, in addition 
to the basic model assumptions, two purely phenomenological ingredients have 
been introduced in order to get agreement with experiment: i) energy s was 
used in the scaled (s/sq) form (with sq being free parameter, this works the 
same way as inelasticity) and ii) the so-called ’’residual function” factor * was 
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postulated {t = —{pi —Pj)^ and (3 being a free parameter) to cut the transverse 
part of the phase space. Therefore sq and (3 were the only relevant parameters, 
e. For our purpose this is sufficient and there is no need to use more sophisticated 
approach exploring escort probabilities formalism, see [ca (for most recent 
discussion of different constraints and their meaning see [na and references 
therein). 

/. Strictly speaking in [na it was shown only for fluctuations of 1/A given by 
gamma distribution. However, it was soon after generalized to other form of 
fluctuations and the word superstatistics has been coined to describe this new 
phenomenon [HU- Another generalization of this idea can be found in [El- 

g. Notice that the y space is limited to y S {—Ym^Ym) where (with pT = 
yjm? + {pt)^ being the mean transverse mass and M' = M — {N — 2)pt ac¬ 
cessible kinetic energy) the limits are given by 


Ym = In {M'/i2pT) [l + (1 - } 


h. It was jus a coincidence that at ISR energies this condition was satisfied. But 
because such a behavior of N as function of energy was only some transient 
phenomenon, there will never be scaling of this type at higher energies, notwith¬ 
standing all opinions to the contrary heard from time to time. 

i. In [E! we have fitted data on pp with q < 1 being the only parameter and 
M = s/s. The q < 1 was cutting off the available phase space playing effectively 
the role of inelasticity K. This result has provided the cosmic ray physicist 
community justification of the empirical formula they used, namely that f{x) oc 
(1 — a • a;)", where x is Feynman variable, x = 2Ej and where a and n are 
free parameters. They turned out to be both given by the parameter q only [ 

[TH]. 

j. In what follows only interactions resulting in the production of particles in the 
central region of reaction are of interest to us, elastic and diffractive dissociation 
collisions will not be considered. 

k. Although this regularity has been observed only for charged secondaries, we 
shall assume here its validity both for the total number of produced particles 
as well as for the fixed number of collisions v. 

l. Notice that such procedure is possible only under the assumption of independent 
collisions. 

771. We would like to stress at this point that our approach to pA using IT concepts 
and tube model differs from that of [EH! because we have only one parameter, 
inelasticity K = R, with both the normalization, shift of the momenta and 
’’partition temperature” being fixed by it, whereas in [ | 2 H| they are all free 
parameters. 

n. We shall not pursued further this problem, which in our opinion can be investi¬ 
gated in the spirit of IT only when the same experiment will provide data both 
for different and well defined values of v and for (z/), i.e. averaged over different 
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V (such possibility was apparently under consideration by NA49 Collaboration 
[ED). We are also aware of potentially interesting new data from RHIC on 
dAu collisions [ |22 and recent attempts of their description [isa and we plan 
to address this issue using IT approach elsewhere. 
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